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Abstract 

The propagator of a scalar field on a stationary slowly varying in space gravitational background is 
derived retaining only the second derivatives of the metric. The corresponding one-loop effective action 
is constructed. The propagator and the effective action turn out to depend nontrivially on the Killing 
vector defining the vacuum state and the Hamiltonian of a scalar field. The Hawking particle production is 
described in the quasiclassical approximation and the quasiclassical formula for the Hawking temperature is 
derived. The behaviour of the Unruh detector on a curved background is considered and the quasiclassical 
formula for the Unruh acceleration determining the Unruh temperature is derived. The radiation reaction 

,— ^ problem on a curved background is discussed in view of the new approximate expression for the propagator. 

O ..... 

The correction to the mass squared of a scalar particle on a stationary gravitational background is obtained. 

This correction is analogous to the quantum correction to the particle mass in a strong electromagnetic field. 

For a vacuum solution to the Einstein equations, it is equal to minus one-fourth of the free fall acceleration 

squared. 
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I. INTRODUCTION 



In calculating the the propagator of quantum fields on a fixed curved background, we encounter 
severe technical problems (see for review, [T]) such that the exact expression can be obtained only for 
a very limited class of spacetimes possessing a large symmetry group. However, for many physical 
problems we need an approximate expression, at least, under the assumption that the gravitational 
field varies slowly. This is the so-called low energy approximation [HE]- In the present paper we 
shall obtain the explicit expression for the propagator of a scalar field on a stationary slowly varying 
in space gravitational background retaining only the second derivatives of the metric. 

The essence of the problem in constructing the approximate expression for the propagator in 
the external gravitational field is related to the fact that the background field stays at the second 
derivatives of the Klein-Gordon operator. In other words, this field enters the principal symbol of 
the operator. So, in the leading approximation when the quadratic part of the symbol is taken 
into account, the metric field should be considered as a constant matrix. Of course, such an 
approximation is inappropriate for sufficiently large point separation of the propagator or for the 
problems where the next to leading corrections to the propagator depending on the derivatives 
of the metric are relevant. The Hawking particle production |3|, the Unruh effect @], and the 
radiation reaction problem [5HS] on a curved background are among such physical problems. We 
shall derive the next to leading correction to the propagator employing appreciably the stationarity 
of the background metric. 

It was shown in [9] (see also |10| HT] ) that for stationary backgrounds the Killing vector 
appears in the effective action in a nontrivial way, that is the expressions involving it cannot be 
rewritten in terms of the local expressions involving the metric and the curvature alone. This fact 
can be used to measure experimentally the variation of the Killing vector (up to multiplication by 
a constant) from point to point. The present paper continues the research in this direction. We 
are about to investigate in detail how the Killing vector field influences (determines) the dynamics 
of the matter fields. The first object for such an investigation is, of course, the propagator of 
a quantum field. Making use of the explicit expression for the propagator, we shall derive the 
analog of the Heisenberg-Euler |12[ IT3] effective action which, as we shall see, depends nontrivially 
on the vector field We shall also obtain the quasiclassical formulas for the Hawking particle 
production [3J and for the Unruh temperature |1] of a detector moving in a stationary gravitational 
field. The quantum correction to the mass of a particle will be also found. This correction is not 
due to the Higgs mechanism, but a direct one, and analogous to the correction to the mass of a 
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charged particle in a strong electromagnetic field |14| . On the vacuum solutions to the Einstein 
equations {R^ v = 0), the correction to the mass squared is found to be equal to minus one-fourth 
of the free fall acceleration squared. The radiation reaction effect on a curved background as a 
mean to measure the direction of the Killing vector will be also discussed. Of course, one may 
contrive less sophisticated methods to measure the various contractions of the Killing vector |1 51 — 
I18j . We address just these problems as their solutions follow immediately from the form of the 
propagator. The mention should be made that the above two problems concerning the Unruh effect 
and the radiation reaction effect are widely discussed in the literature arriving at the contradictory 
conclusions (see, e.g., |19| I2U]), One of the aim of the present paper is to provide a more solid 
mathematical background for these considerations on a curved background. 

The paper is organized as follows. In Sec. [TTJ the derivative expansion of the Feynman propaga- 
tor is considered near its diagonal up to the second derivatives of the stationary metric field. The 
stationarity of the metric is essentially employed and the generalized Schwinger-DeWitt technique 
|21j adapted to the (3 + l)-decomposition is used. We shall also see that the simple derivative 
expansion has a large degree of ambiguity and cannot be used to restore the propagator unambigu- 



ously. In Sec. Ill the derivative expansion is resummed making use of the Schwinger approach [13j, 
and the low energy expansion [U |2] of the propagator is constructed thereby. Here we shall derive 
the explicit expression for the propagator (positive-frequency function) as the double integral over 
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the frequencies of the modes and an auxiliary variable - the Fock proper-time [22] • Section IV 
devoted to a further simplification of the expression for the propagator under the reasonable addi- 
tional assumptions on the fields and the point separation. Several physical implications following 
immediately from the form of the propagator will be discussed. We shall obtain the spectral repre- 
sentation of the propagator and its representation in terms of the proper-time. The expression for 
the one-loop correction to the effective action will be also derived. In conclusion, we shall discuss 
some implications of the results of the present paper. In appendix |Aj several formulas concerning 
the calculus on the spacetime with the Killing vector are presented. There, in particular, the ex- 
pressions for the 4-vectors of the gravitational and inertial forces constructed with the help of the 
Killing vector field £ M will be given. Some useful expansions needed in the course of the derivation 
of the approximate expression for the propagator are collected in appendix [B] 
We shall use the conventions adopted in |21| [23] 

for the curvatures and other structures appearing in the heat kernel expansion. The square and 



round brackets at a pair of indices denote antisymmetrization and symmetrization without 1/2, 
respectively. The Greek indices are raised and lowered by the metric g^ v which has the signature 
+2. Also we assume that the metric possesses the timelike Killing vector £ M 



c&w = o, e = g^ee < o, 



(2) 



that allows us to make the decomposition \2 



ds 2 = g^dx^dx" =: ^(g^dx^) 2 + g^dx^dx", 



(3) 



where = ^/^ 2 is a one-form dual to the Killing vector (the Tolman temperature one-form). In 
the system of coordinates, where £^ = (1,0,0,0), we have the relations 



9ik9 kj = S(, f det^ = g, g 00 - m,g 0l g 0j = (goo)' 1 , 
o o 1 



9iw 



gij 



9i09j0 
900 . 



g 00 - goo 1 9 0J 



g 



i0 



g 



(4) 



The the Latin indices corresponding to the space are raised and lowered by the positive-definite 
metric g^. The curvatures associated with this metric will be distinguished by the overbars, e.g., 
R. The system of units is chosen such that c = H = 1. 



II. ASYMPTOTIC EXPANSION OF THE PROPAGATOR NEAR THE DIAGONAL 

Our aim in this section is to obtain the derivative expansion of the Feynman propagator near its 
diagonal up to the second derivatives of the stationary metric field g^ u . We shall take exactly into 
account the symmetry of the background metric with respect to translations in time and employ 
the so-called generalized Schwinger-DeWitt technique |21| adapted to the (3 + l)-decomposition. 
Hereinafter we shall work in the reference frame where the Killing vector is straighten = 
(1,0,0,0)). Besides all the formulas presented in this section are valid for an arbitrary spacetime 
dimension (d + 1). 

Let 



H(x,y) = (V 2 x -m 2 



5(x - y) 



\g\V*(x)\g\V*(y) 

= \g\~ l/ \x) 



- 1 ^(x)d^\g^d u \g\- 1 /\x) - m 2 



S(x — y) 



(5) 
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be the bi-scalar kernel of the Klein-Gordon operator. Supposing the metric does not depend on 
time t, we can write 



H(x,y) =: |C 2 r 1/4 (x) / %-^H(co; x, y) l^ 2 1 " 1/4 (y ) , 

2 1 1 1 X( _ ^ W 

" +( y i + iug . ) 2_ rn 2_^%- 1 A h i h i d[X yj 
5 z 4 



» 1 /4( x ) g V4(y)' 

where x = (i,x) and y = (0,y), the Levi-Civita connection Vj is constructed in terms of the metric 
(jij, and we have introduced the notation ha '■= In v/|£ 2 |. Notice that the one-form gi plays the 
role of a £7(1) gauge field with the charge u) (the energy). The kernel H(co; x, y) is a bi-scalar with 
respect to the spatial diffeomorphisms. Then the positive-frequency function defined with respect 
to the standard vacuum state for the stationary backgrounds (in black hole physics it is known as 
the Boulware vacuum [25J) is written as 



OO 



D^\x,y) := -mx)ij>{v)) = "* / da;e-- t -|e 2 r 1/4 (x)(x| ( J(ff(a;))|y)|e 2 r 1 / 4 (y), 

j° (7) 
{K\S(H(u))\y) = J g(x|e^)|y> =: J gG(u>, s; x, y), 

where i_ := t — iO and the integration contour in the s-plane passes a little bit lower than the real 
axis. The Feynman propagator is expressed in terms of the positive-frequency function: at t > 
the propagator coincides with it, and at t < we should make the replacement, 

t-^-t + iO, xoy, (8) 

in D^ + \x,y). So, it is sufficient to calculate the positive-frequency function only. 

To this end, we can apply the generalized Schwinger-DeWitt technique to the heat kernel 
G(ui, s; x, y) assuming that the points x and y can be connected by a unique geodesic of the 
metric Making the derivative expansion, it is convenient to use the covariant gp-symbol of the 
heat kernel 1 26 1 



G = ^is(-^+D 2 -m 2 --V i h l - l h t h^ J _y^ e ^'a (x,y) = 

d d k' /— itf' a^-is^+k^ +m 2 + \V i h i + \h i h?) , > 

g'e i »* / 2 4 ao(x, y)x 



(27r) d 

, s is(£+k^ a^a* + ±V i h i +ih i h i ) ^(-^+(A+^'a l , l ) 2 -|V'? ll -ifc I / l l ) , , fn . 

xao(y,x)e « J 2 4 e « ! 4 ao(x, y), (9) 

where the primes denote the quantities and derivatives referring to the point y, and ao(x, y) is the 
geodetic parallel displacement operator with respect to the connection D{ = Vj + iuigi, 

ao(x J y) = e- iw ^ <fai *. (10) 



The integral is taken along the geodesic of the metric g~ij connecting the points x and y. 

The derivative expansion is obtained expanding the last line in ^ in a covariant Taylor series 
near the point x. At first, all the exponents in this line are expanded in the series. After that the 
covariant derivatives D 2 are carried through the expression to the right to act on ao, whereas the 
derivatives of ao are expanded in the covariant Taylor series by the formulas (4.26)-(4.29) in |21| . 
In the notation of |21j 

11^ -> iuifij, (11) 

where := du,9j\- The first exponent in the last line of Eq. (|9| provides considerable cancellations 
in the expansion - all the terms, which do not contain the derivatives of the expression staying 
in the exponent, are gone away. The exponent in the second line effectively sums such terms (see 
for the proof [27, 28J). As a result, we obtain the series containing the nonnegative integer powers 
of k 1 ' . Then the Gaussian integral over k 1 ' is performed. In fact, one can replace k 1 ' by —id/ddii 
acting on the result of the Gaussian integration of the expression in the second line in ^ 

^^ e ^ 1 ^>V^(m' + g + ^ + jM«)^ ( ^ (12) 
(4:iris) d ' 2 

After rather lengthy but straightforward calculations one derives 



^^e--- (m2+ # + ^ l ^^^\ (x,y) jl + %o*a> + %R- 
(Ams) a l z [4 6 

V l hi - -hih 1 ) - —V J fjia l - —fijf J + —-^hio* 



is 



s 

3 



V3 1 3 / 3s JJl 2 JlJJ ' s 

-.A . ,3 o, ,2 , 1 s , ,2 2 -, 2 4 ^ 

(13) 

up to the second derivatives in the background fields, all the fields being taken at the point x. 
Notice that in this approach the order of the term in derivatives is counted as the total number of 
derivatives entering it. For example, the term 

(hih*) 2 (14) 

is of the fourth order in derivatives. 

The heat kernel G(ui, s; x, y) is a unitary operator with respect to the standard measure g 1 / 2 or, 
equivalently, it is Hermitian with respect to this measure when s is purely imaginary. It is desirable 



to preserve this symmetry in the approximate expression for the heat kernel. The expression (13) is 



not Hermitian at s = ir since the points x and y enter asymmetrically to it. Therefore, we expand 



(13) in the vicinity of the point p(x, y) lying at the middle of the geodesic connecting the points 



x and y (the midpoint prescription) and retain only the terms which are no more than the second 



order in derivatives of the fields. An inspection of the expansion (13) shows up that we only need 
to expand £ -2 in the exponent and the last term in the second line, while the fields entering the 
other terms are simply taken at the midpoint. This gives 



_ q (x,y) is-i s (m 2 +^ + iV'/i,+i/ i ,/ l ') 
(AiTis) d / 2C 

+ ^R+ii^r^h _ 2/a 7 w - - 

6 12 £ 2 v 3 31 6 



1 + —R ij a i a j + 



e 

is 
~3 



2V l hi 



Ah t K l 

,3 



V'./>' 



f. . m 

2 J % 3J 



3 r W 4 



h (15) 



where we also expand the bi-scalar van Vleck determinant A(x, y) (see Appendix |BJ) . 

Now we should substitute the expansion obtained to the integral ([T]) and integrate it over u and 
s. However, at this point we have a large degree of ambiguity how to represent the substituting 
expression, and this ambiguity can change substantially the resulting positive-frequency function. 



For example, we may expand the exponent in ( 15 ) and neglect the higher order terms or, conversely, 
collect some terms in the preexponential factor to the exponent similarly to what we did with the 
term w 2 /£ 2 . The naive derivative expansion, which we carried out in this section, does not tell us 
what the right form of the integrand is. We shall address this problem in the next section, while 



here we establish the relation of the expansion (15) with the standard one |21[ 123] written in terms 
of the geodetic interval of the spacetime rather than a. In such a way, we shall make a nontrivial 



check of the correctness of the expression (15). 

In order to obtain the standard asymptotic expansion of the Feynman propagator, it is convenient 
to employ the Schwinger representation for it 

ro ° doj 



D(x,y) := -i(T{<Kx)<t>(y)}) 



2?r 



iut I £2 1-1/4 



(xKxKi^ + iorV^r 174 ^), 



(16) 



lm 2 — >m 2 — iOi 



poo poo 

(x|(i/(a;) + iO)- 1 |y) = / ids(x|e^ (a;)+i0) |y) = / idsG(uj, s; x, y) 

Jo Jo 

instead of using Eq. 0. Then we can formally integrate over uj the expansion (PL5J) as the integral 
of a Gaussian type. This leads to 



\er 1/4 G\$ 



/2i — 1/4 



(4vris)( d + 1 )/ 2 



% 

48s 



£ 2 T 2 (Vify - 2hihj + ffikf k j )a i (t j + 



4^4 . 2 



+ 1 + 



l r 

12 



FL 



hihj 



2rj-\2 1,2 



+ 



IS 



R - 2h 2 - 2V l h; 



4 3 



(17) 
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where T := t + /* dx l gi is a scalar with respect to the general coordinate transformations 
Bearing in mind that = at fi = in the adapted system of coordinates (i.e. £^cx M = 0) and 
applying the formulas given in Appendix [X] we can write to the same accuracy 



|£2|-1/4 G |£/2|-1 



/4 



48s 



(47ris) (d+1)/2 I 

+ 1 + ^R^ire + ^)(rr + o v ) + jR\ 



exp 



^-{?T 2 + 2a- ^[£ 2 (T^ + f^a u f - (V^K - 2h^h v )^a u ]} - ism' 



+ ^ATe + ^){n u + ° v ) + jr}. (is) 



( 47ris )(d+l)/2 

Up to the terms of the second order in derivatives, the quantity staying in the curly brackets in the 
exponent is the geodetic interval squared, and 



(19) 



to the leading order in derivatives. 
Indeed, it is sufficient to consider 



X 2 (x,y) :=f(p)T 2 + 2a 



on the arbitrary worldline x^(r), where r is the natural parameter 

^(x(T))^(r)x"(r) = -1. 



(20) 



(21) 



For definiteness, we can take x M (0) = y = with zero being the origin of the Riemann normal 
coordinates of the metric g^j (see for details |29|). To fix the frame in the spacetime uniquely, it is 
convenient to use the Fock gauge |24| [26) [30] 



9i (x) x i = O gi {x) = 7^iy xjl • • • xj " d n • ■ ■ d j n -Jjm = ^ Sji + ^u^V h f hi + . 



n=l 



(22) 



where x l are the coordinates of the Riemann frame. This gauge is equivalent to 

dx^i = 0, (23) 
for any x l of the Riemann normal coordinates, where the integral is taken along a straight line, and 



so T = t. Note that the higher terms of the expansion (22) depend on the curvature [29], but they 
will be irrelevant in our case as they are of the higher order in derivatives. Then expanding the 



expression (20) in r and making use of the relations (21), (22), one arrives at (see Appendix [B 

o 1 



X z 



12 



(v uU y T=0 T* + ... 



(24) 



up to the terms of the second order in derivatives of the metric g^. The Killing vector disappears 
from the expression. In a free fall, X 2 = — r 2 and, consequently, with the accuracy we work, X 2 



coincides with the geodetic interval squared. Thus we see that the expansion (18) agrees with the 



standard asymptotic expansion of the heat kernel near the diagonal |21| I23| . So, the expansions 



for the propagator following from (16) are also the same up to the terms of the second order in 
derivatives. 

III. RESUMMATION OF THE EXPANSION 

In the previous section we saw that the resulting expression for the positive-frequency function 
-t)( + ) (x, y) depends severely on the rearrangements of the derivative expansion. The very fact that 



we was able to integrate the expansion ( 15 ) over the frequency uj is a consequence of the resummation 
of the derivative expansion, which we have done raising the potential term w 2 /£ 2 to the exponent. 
Therefore we need a more reliable procedure to obtain an approximate but adequate expression for 
the propagator. This procedure should sum an infinite number of terms of the derivative expansion, 
say, all the terms containing the fields at the given point, their first, second derivatives and not the 
third derivatives and higher. This is the so-called low energy expansion of the heat kernel [H|2]. For 



example, in this case the term of the form (14) must not be neglected. In this section, we obtain 
the leading term in the low energy expansion of the heat kernel slightly generalizing the standard 
procedure |13] to the curved spacetime. Now we put d = 3. 

To begin with, we consider a simple model in a flat Euclidean space with the action 

S[x(t)\ = J^dr^x 2 - Ai(x)x l + itfysV), (25) 

where Ai = x J fji/2, the field strength matrix is constant and skewsymmetric, Eij is a constant 
symmetric matrix. This is a general quadratic model and later on we shall see how to use it to 
construct a reliable positive-frequency function of the scalar field on a curved background. Further, 



we assume 



[E, f] = / = -tHv^ =: -iHaij, E = Xivft^ + X 2 vfvj =: X lSij + X 2 vfv 2 . (26) 

The vectors vj, vh and vf are orthonormal with respect to the standard Hermitian scalar product, 
the overbar denotes complex conjugation, the vector v 2 having real components. In this case, the 
equations of motion for this model can be readily integrated 

x = [{Qie- iu>+T +C2e- iul ~ T )v l +c.c] +(c 1 e^ 2T +c 2 e-^ T )v 2 , u± := ^(H± sj H 2 - 4A X ), (27) 



where £1 and £2 are the complex numbers, c\ and C2 are real, and we assume for definiteness that 
A2 > and Ai < 0. This is indeed the case for the gravitational field (see below). The constants of 
integration are uniquely determined by the boundary conditions a;(0) = y and x(s) = x. Then the 
Hamilton-Jacobi action takes the form 

S = ^{xiis^is) - xi(0)i*(0)) = \{y l s ijy i + x's^WH 2 - 4A X ctg | ^H* - 4A X - 



\JH 2 - 4Ai . , s TT . , i . s TT . 

--(x Sijy J cos -H + ix a,ijy J sin -H)+ 



2sinf v / ^ 2 -4Ai 



cthsVA 2 [KV) 2 + KV) 21 



^x^V = -y^/-P-±E C tg y-p-AEy+ 
shsyA2 4 2 



+ 7^V / -/ 2 -4gctg-y / -/ 2 -4gx-j/ ^ /2 ; ^' La;. (28) 
4 V S 2 V 2sinfV-/ 2 -4f? 



Shifting the variables in the initial action (25), it is easy to generalize the above result to the case 



where the Lagrangian contains the additional terms Eq and b{X l . And so, for the system with the 
Hamilton function 



H(p, x) = {pi + Ai) 2 -E Q - b i x i - ^ /•:,,. 



(29) 



we have the Hamilton-Jacobi action 



S = \{x - y)xctg sx{x -y) + y^(ctg sx - -?—)x - S -bE~ l b + \bE- 1 f{x - y) + sE , (30) 
4 2 sins/ 2 2 

where x := \/—f 2 — 2E, while Xi = x% + E^bj and jji = y% + E~^bj. Note that the factor 1/2 



is absent at the kinetic term in the Hamilton function (29). Therefore we have to stretch the 



proper-time and the potential in (25) accordingly in order to get (29) and (30). Also notice that 



the Hamilton-Jacobi action (30) contains the inverse matrix J5 f • , but it has a finite limit when E, 



'j 



becomes degenerate as seen from (30) expanded in s 
{x-y) 2 1 

„3 r 



s 



r l l l 

+ 7,xfy + s —{x- y)f 2 (x - y) + -(xEx + yEy + xEy) + -b(x + y) + E 
4s 2 L 12 o 2 



— xEfy — — 
6 Jy 12 



1 

15 



(x-y)(f 2 + 2E) 2 (x-y) + xE 2 y +... (31) 



All the terms with E-- 1 cancel out. The van Vleck determinant is written as 



det 



d 2 s 



-r = (-2s)- d det 
<j 



, , sne s f , „ . j , /sinsxx- 1 



(-2s)- d det - 



1 sin sx 

- m 

2 



where we have used the unimodular property of the matrix e s f . 



(32) 



10 



Now we return to our problem. The leading (Gaussian) contribution to the low energy expansion 
of the heat kernel is obtained |13j if we expand the Hamiltonian ^ in the momentum and coordinate 
operators up to the second order and throw away the higher terms. Then we solve the Heisenberg 
equations for the retained quadratic Hamiltonian and find the evolution operator (the heat kernel) 
in the Fock proper-time |13[ l22l 131] , As well-known, this problem is essentially classical for the 
quadratic Hamiltonians, and the quasiclassical answer, 



G( U ,8;x,y) = <x|e- is M^%) = [(-2™)"'* det ^-j ' (33) 

for the Green function is the exact one. Here tilde reminds us that the Green function is the 
bi-density rather than the bi-scalar as in Eq. Q. Formula ([33} can also be obtained from the 



d 2 S n 1/2 , 



path-integral representation of the heat kernel after a Gaussian integration over the fields x l {r) 
(see, e.g., [2311231321). 

The Hamiltonian for our system is 

-H = g- l, \pi + 9i)V§9 ij (Pj + 9j)r 1/4 + \Vhi + \hW + ^ + m 2 = 

= {Pi + 9i)9 ij (Pj + 9j) + \di{g ij dj In + i$ In y/gfldj In y/jj + + ifctf + ^+m 2 . 

(34) 

Of course, we want to derive such an approximate expression for the heat kernel associated with 
this Hamiltonian that preserve all the symmetries of the exact evolution operator: it should be the 
kernel of a unitary operator, it should be the bi-density function which is invariant with respect 
to the gradient transformations of the field g^. The last property guaranties that the correspond- 
ing positive-frequency function will be generally covariant under the spacetime transformations. 
Besides, we want that already the leading (Gaussian) approximation give us the most exact ap- 
proximation that we can achieve for the heat kernel in the case of the slowly varying fields g lJ and 
gi. So as we need not to evaluate the higher order corrections to it using the perturbation theory, 
or reduce this work to a minimum. 

To this aim, we adopt the following strategy. For any given points x and y connected by the 
geodesic, we pass to the Riemann normal coordinates with the origin at the midpoint p(x, y). 



In this frame, we expand the functions entering the Hamiltonian (34) in a Taylor series. This 
automatically gives us the covariant expressions at any finite order of the expansion. In particular, 
using the formulas presented in |27[ [29] and in Appendix p] we find 



^dttfldj In ^) + \di In hx^/§ = ~R- ^ViRx' + ^xV + . . . 

t . \ ( _ E> fjk __ _ r>kl f=> _ mnk f> _ fb ^vv f>\ 

ij ■ — ^ \ r^ rl ik rl j g-^ ^kilj g * ^jmnk ^ v Jr H] ^ Vijlh}. 
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(35) 



In passing from the bi-density to the bi-scalar, one should multiply the former by 



F 1/4 (*)r 1/4 (y) = A 1/! W), 



(36) 



where we have used the relation |27] 

A(x,0) = A(0,x)=r 1/2 (x), 



(37) 



and the composition property of the covariant van Vleck determinant. Also we need the expansion 



u} 2 C (x) = u C 2 (l - 2h iX l + (2hihj - Vih^x'x 3 + ...)=: w 2 (<f 



W-*-iEgW). (38) 



b u 'x 



Hereinafter, all the fields and their derivatives are assumed to be taken at the point p, unless 
otherwise stated. As long as we use the midpoint prescription, a 1 = 2x l = —2y l . 



Now, in the Hamiltonian (34), we retain the terms which are at most quadratic in x and p 



(keeping the order that they are written in (34)) and obtain the Gaussian approximation for the 
heat kernel 



G(u),s;x,y) 



M 



AV2 (xy) /gjngrpi 

m 2 - -R+-V i h i + \h 2 
6 2 4 



1/2 



iSo(s,ai)—isM 2 



m 2 -l^R + ^R^^-lh 2 , 
e sulf s ,a 



S : = ]axctg(sx)a bE'^x^tgsx - ^—)(~ + E~ l b) - ~-bE l />+ 

4 2 2 smsx 2 2 



(39) 



+ ^bE- 1 fa-s^, 

where for brevity we omit the matrix indices, a = a% in what follows, x is defined as above with 
the replacement / — > uf, and 



E. 



, ,2 p(0) , p(2) 
U E ij + E ij ' 



u 2 b^+b?\ 



b^ = V, 



l -{ R+ eR^9 v ) + \h 2 



(40) 



Expanding the expression for G in the number of derivatives as it was done in the previous section, 



one can convince oneself, using formulas (31), (32), that all the terms of the expansion (15) are 
reproduced save the term 



(41) 



It seems this term cannot be obtained from the Gaussian approximation. Such a term may come 
from the divergence of the gauge field gi that appears when one uses the gp-ordering of the Hamil- 



tonian (34). However, the gp-ordered Hamiltonian truncated to its quadratic part is not Hermitian 
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and the evolution operator is not unitary - the property which we desire to preserve. Thus, substi- 
tuting the expression obtained to Q, we arrive at the central result of the present paper 

' ' 2tt |^ 2 | 1 /4(x)(47ris)' i / 2 |^ 2 | 1 /4( y ) V sx J 







(42) 



where t = T or the Fock gauge (22) based at the midpoint is implied. The higher derivative term 
(41) could be "exponentiated", but we shall not investigate this possibility. 

Some remarks are in order. Due to the noncommutativity of x and p in the kinetic part of the 



Hamiltonian (34), there is an ambiguity how to single out the quadratic part from it. For example, 



one could use the Weyl ordering, where the first term in (34) is rewritten (without the gauge fields) 



Pi9 ij Pj = ' (PiPj9 ij + 2Pt9 ij Pj + 9 ij PiPj) + \dijg ij - (43) 

Then one extracts the quadratic part from the Hamiltonian and obtains the same associated classical 
equations of motion as above, but with the additional correction to the potential. In particular, 



instead of —R/6 in (35) one will have 



ifl+l^i(0) = -ifl. (Hi 



It is this correction to the potential term which is argued by DeWitt as the correct one (see [23]), 
when one uses the path-integral representation of the heat kernel with the midpoint prescription. 
However, with this potential term, one needs to perform the two-loop calculations in order to 
reproduce the well-known asymptotic of the heat kernel on the diagonal even at the first power of 
s. The potential with such a property is, of course, quite unsuitable for our purpose outlined above. 
In fact, changing the ordering prescriptions, one can obtain any number at the scalar curvature in 
the potential. The conditions, that we impose on the approximate heat kernel, fix the prescription 



unambiguously. It was proven in |27| that the correction of the form ( 35 ) to the potential term 
sums all the terms of the derivative expansion of the heat kernel, which contain R at the given 
point. 



In order to provide a tighter connection of the positive-frequency function ( 42 ) with reality, let 
us write the quantities entering it in the weak field limit for the vacuum solutions of the Einstein 
equations (for details, see j2j and Appendix |A| 

1,(0) ~ r 9„ l(2)_ r 9 n p(0) _ r g 

t ^ ^ 3/*^ 39 

2r„ 3 2 3r 2 

fij ~ ^k( Mi i + 2 n [i M j\k n k), [E, f] ~ w —^n^M j)k n k , 
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where m is a total mass of the gravitating object, Mij is its angular momentum, r g is the 
Schwarzschild radius, r is a distance from the gravitating object, and n{ = x%jr. In the sec- 
ond line, the Schwarzschild metric was used for the calculations. We see that the tensor E^ has one 



positive and two negative eigenvalues as in (26). The relation [E, f] = 0, which we also assumed, 
is fulfilled when the vector of the angular momentum is parallel to the vector of the gravielectric 
force (see Appendix [A] for the definition) . The above calculations can be easily generalized to the 



case when the tensor /j,- has a canonical form (26) and the tensor Eij is diagonal in the same basis. 
One should set Ai 6 C in the solution (27), and the matrix e s f in the Hamilton- Jacobi action (30) 
is to be placed near x. Nevertheless, further we restrict ourself to the case [E, f] = 0. 



IV. PHYSICAL IMPLICATIONS 



In this section we shall discuss some physical implications of the result obtained above. The 



crucial point is, of course, the fact that the positive-frequency function (42) depends nontrivially 
on the Killing vector. This dependence stems from the dependence of the Hamiltonian of the scalar 
field on the Killing vector of the gravitational background. The positive-frequency function is an 
observable and so it is just a matter of the physical erudition to extract the dependence on the 
Killing vector from it. We shall find the various expressions for the propagator of the scalar field 
which are valid in the different regimes: at sufficiently large and extremely small point separations. 
Besides, we shall derive the exponentially suppressed contributions to the one-loop effective action 
which depend nontrivially on the vector field This action is induced by the polarization of 
the vacuum of scalar particles and represents the analog of the Heisenberg-Euler action |12|. [T3] in 
quantum electrodynamics. In particular, we shall obtain the quasiclassical formulas for the Hawking 
particle production [3] analogous to Schwinger's formulas for the particle production in the constant 
electromagnetic field. The Unruh effect and the radiation reaction problem on a curved background 
will be also addressed and the formula for the acceleration determining the Unruh effect will be 
derived. 



The structures appearing in the heat kernel (39) fall naturally into the pieces denoted by the 



indices and 2. Let us estimate their ratio. Introducing the characteristic scales in the weak field 
limit 

9 i ~r" 1 ~L" 1 , hi^e/L, Rij ~ e/L 2 , etc. (46) 
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where e = (1 — £ 2 ) ~ r g/L <C 1, we see that the inequalities 



u 2 m»\b^\ 



are equivalent to 



lo 2 L 2 > 1. 



Near the horizon, e«l, the estimations (46) changes, 



etc. 



and the conditions (47) become 



eVx 2 » i, 



(47) 



(48) 



(49) 



(50) 



which is valid in the weak field limit as well. Of course, certain combinations of the fields <?j, /y, 
R % jkl, and their derivatives may be unexpectedly small or vanish on the particular solutions to the 
Einstein equations, but their orders of magnitude are not larger than given above. Loosely speaking, 



the inequality (50) says that we consider a propagation of the wave packet consisting of the modes 



with the wavelengths much smaller than the distance to the gravitating object. Therefore, in the 
range of applicability of the approximation of slowly varying fields, which we assume from the 



outset, the inequality (50) has to be hold and the estimations (47) are fulfilled. 



Now we evaluate the different types of contributions to the integral ( 42 ) . Stretching the proper- 
time s — > s/co and redefining 



E: 



(51) 



we can write the integral in the form 

ds AV2( X} y)(i + a vVji<r76) 



L>(+)(x,y) = -i 



2vr|e 2 | 1 /4( x )(47ris) d / 2 |e 2 | 1 / 4 (y) 

dww d/2-1 det[ 



/sinsx\ - 1 / 2 



o 



SK 



jLL}(So(s,a-i)-t-)-isM 2 /oj 



(52) 



where So is defined as in Eq. (39) with uj = 1. In the representation (26), 

-J2 



S = ^rVH 2 - 2Aifctg sVH 2 - 2Ai + 



cos sH 



+ ^ VH 2 - 2Ai (ctg s^jm - 2Ai 



sinsV-tf 2 - 2Ai 
cos sH 



at 



+ 



7 2A 2 th-V2A 2 + 



sin sy/ H 2 
sin sH 



-z-(a x b y - a y b x )Uii 2 -2Ai-. 

2Ai \ smsV-n — 2Ai 



H 



+ -f2V2A 2 cth 

s ,b 2 "|| _s_ 
, e2 



2A 2 



(53) 



2 v Ai A 2 ' 



where t; 2 is assumed to be directed along the z-axis and the projections to this axis are denoted as 
parallel. 
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A. The ^-representation 



At first, we take the integral over s. The integrand has the three types of singular points in the 
s-plane: (a) the essentially singular points on the imaginary axis s = z7rn/-v/2A 2 , (b) the essentially 
singular points on the real axis s = im/y/H 2 — 2Ai, and (c) the branching point at the origin. The 
points of the type (a) are responsible for the Hawking particle production, the points of the type 
(b) describes the vacuum polarization effects, and the point (c) gives the major contribution to the 
propagator. None of the contributions from these points can be evaluated exactly. However we can 



make it approximately under the assumption that the estimation (50) holds. We shall expand Sq 
near the singular points retaining only the leading terms of the Laurent series. The higher terms of 
this series give a negligible contribution. This is easy to see if one makes the variable s dimensionless 
stretching it once more as s — > s/u). 
Introducing the notation 

V# 2 - 2Ai , 



12 : = 



/2A 2 
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H 



/2A 2 ' V# 2 - 2Ai ' " m ' = VfP 

the expansions near the points (a) can be cast into the form 



2Ai 



Ih 



2 • = 



H 

7IV 



(54) 



So 







~A * 

Ax 


. 8 


1 t- 




2A^ 


by ~ 



^ dnrnl H 2 
sh 7rn/i2 
-sh -Knlu2 



sh 7rn/i2 



H + 



2A 2 



+ ^rky/H 2 -2Aifcth7rnZ 



'12 



ch imlH2 
sh-7rnZi2 



1. 



(55) 



where x tends to zero and n is an odd number, and 



So 



X 2 x 



\/H 2 - 2Ai(cthvrn/i2 + 



ch7rnZ//2 
sh iml\2 
sh imlH2 



+ -^V^ 2 -2Ai(cthvrn/i2 
2A-^ V 



ch %tiIh2 
sh ixnlyi 



H + 



2A 2 (r + 2^ 



(56) 



when n is an even number. The mass term and the preexponential factor become in these cases 
mnM 2 



sM l 



_ 3/2 /sinsx\-i/2 _ e -™(n+i)/2^jp _ 2Ai(2A 2 ) 1/4 

S CIGTj I 



(57) 



./2A2 ' "\ sx J shvrn/i2sh 1 / 2 (a; v / 2A2") 

Recall that the integration contour in the s-plane goes along the real axis below the singularities 
lying on it. The half-plane (upper or lower), where we should close the contour, can be determined 
analyzing the asymptotic behavior of the integrand at the large n. In the case when 



^H 2 - 2Ai > uj(C 2 + h>E- x b) + M 2 /uj, 



(58) 
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we can close the contour in the upper half-plane. Otherwise, we have to close the contour in the 
lower half-plane and the singular points of the types (b) and (c) do not contribute to the positive- 
frequency function at such lj's. One may say that such modes do not propagate. 



For the energies oj satisfying the estimation (50) and much larger than the Compton wavelength, 



the inequality (58) is fulfilled. Then the contributions of the singularities (a) are suppressed by the 
Boltzmann-like factor 



<'xp|^-!=(£ 2 + \bE l b)} « exp 



irn 



(r 2 + ^ (0) ^& (0) ) 



2A 



(0) 



(or 



(59) 



where A 2 is the eigenvalue of the matrix E^°\ The coefficient at the energy u at n = 1 can 
be interpreted as the reciprocal temperature of the Hawking radiation. Inasmuch as the one-loop 
contribution of one bosonic mode to the effective action reads as 

r$ = -i J J q dujuj 2 D^\uj,x,x), (60) 

these singularities contribute to the imaginary part of the effective action and suppressed by the 



same exponential factor. The spatial measure in the integral (60) is that measure with respect to 



which the mode functions are orthogonal. To evaluate the integral over s, we use the formula 

dz 4 v / 2^r(5/4) 



h shV^zv^Aa") 



/ 2A7r(3/4) 



(61) 



where H is the Hankel contour that runs from +oo a little bit higher than the real axis, encircles 
the origin, and then goes to +oo a little bit lower than the real axis. With this integral at hand, 
we can write the contributions from the singularities of the type (a) for odd n (n = 1 is the leading 
contribution) in the form 



D(+)( W ,x,x)| (a) 



r(n/2+l/4) (2A 2 ) 1 / 4 r(5/4) l 12 

\/27r|S 2 | 1/2 r(3/4) shirnlu 



X UJ 



^H 2 -2X 1 (cth nnh 



ibE^b- 



(62) 



The imaginary part of the correction to the effective action coming from these singularities is 
positive as it should be. There exists also the exponentially suppressed contribution to the real 
part of the effective action from these terms. As far as the contributions with even n are concerned, 



the nonvanishing term at 1/x in the expansion (56) leads to a complication in evaluating this 



integral and the resulting expression is rather huge. So, we do not write it here. 
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Near the points (b), we have 



So 



r ^(l + (-l)"cos7mZ m ) + ^(l 



H(-l) n smirnl H1 ( 



-l) n cos nnlHi) 
a 2 b 2 



-l) n 



2Ai 



a x b y - Oyb x ) simrnlHi 



-^) + V2A 2 cthvrn/ 2 i( 



l + 2X 2 2 )+ 



H TTTl 

+ 2X x ^ h v ~ ^AX 1 ~ (-l) n cos7rn/ H i) - ^ R2 _ 



(r 2 + 2 bE- l b), 



and 



sM 1 



irnM 2 
VH' 2 - 2Ai 



" 3 / 2 det 



/sin sx 



-1/2 (-1)«(2A 2 ) 1/4 



X Sh 1 / 2 (7T?7.Z 



(63) 



(64) 
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The integrand of (52) is a single- valued function in the neighbourhood of these singular points 



because of the degeneracy of the eigenvalue \/H 2 — 2Ai. As a result, the integral over s going 



around these points can be simply evaluated. Assuming the inequality (58) is satisfied, we arrive 
at 



A 1 /2 (X)V )(l + 



l) n (2A 2 



D (+)/ U} x v )| ~ \ Q^/m-2X 1 , 

(b) ~ |e 2 | 1 / 4 (x)(4^)3/ 2 |e 2 | 1 / 4 (y)sh 1 / 2 (vrn/ 2 i) 



I 1 /* 



(65) 



where Sq is the finite part of the expansion (63) at x — > 0. 

The main contribution to the positive-frequency function comes from the origin of the s-plane. 



From the expansion (31) near this point, we deduce 
-24 



o 



O" 



+ s[±*(2f + E)*-±] 



6 J 12 



[b 2 + ^-a{f + 2E) 2 a - \aE 2 a\ +... (66) 
L 15 4 J 



Rescaling the proper-time s — > s/lo and bearing in mind that Sq is multiplied by u> in the exponent 



(52), we see that the terms of the above expansion at the second power of s and higher are small 



provided the condition (50 ) holds. Besides, there are also the saddle points near (after the rescaling) 

E := -l/i 2 + a(2f 2 + E^)a/24, 



u 2 a 2 



m 



:= M 2 -aE i2) a/2A. (67) 



A{E-m 2 /uj 2 ) 

If these points are situated near the origin they will considerably contribute to the integral. Here we 
are interested in the approximation of the propagator for the point separation much larger than the 
wavelength of a mode (for the extreme case of the infinitely small a and t see the next subsection). 
So, we assume that those extremum points are far from the singular point s = 0. This condition is 
fulfilled when 



leW > i. 
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(68) 



Hence, keeping only the terms at s , s , and s in the exponent of (52 ) and the leading contribution 



from the preexponential factor, we obtain the approximate expression 



D {+ \x,y)\ 



ds 



, ,d/2-l _ 2 

J - Juj(f I -t-+sE)-ism 2 /Lo 



(69) 



|e 2 | 1/4 (x)|e 2 | 1/4 (y) J 2vr 7 (4vris) rf / 2 
The integral over s is reduced to the Bessel function J]y 2 for d = 3 and, consequently, is expressed 
in terms of elementary functions. After a little algebra, we find 



do; 



0(c/£7 - m>~ i0Jt - sin 



(70) 



|e 2 | 1 /4(x)|e 2 |V4(y) 7 47r 2 |a| 
The latter integral is the same as the integral for the positive-frequency function in a flat spacetime 
at rh 2 > with the obvious redefinitions. Applying the formulas from the appendix of [33], we 
come to 

-AVZfoy) fsgn(i) 



D (+) 



(c) 



+ 



|e 2 | 1/4 (x) J E 1 /2|^2|l/4( y ) [ 4vr 
irh9(\) 



-S(X)+ 



[iVi(m\/A) + isgn(t)Ji(mVX)] + im6 } J^ K, (m 



sgn(t) 



5(A) 



+ 



4vr 2 A 8vr 



In 



A) 
mjA) 1 / 2 



(71) 



m 
16tt 



sgn(t)0(A) 



|^ 2 | 1 /4(x)f;l/2|^2|l/4(y) L 4vr 

A = -A >2 (x,y) :=^" 1 t 2 -a 2 , 

where the last approximate expression is the expansion of the positive-frequency function at the 
small mass rh. Of course, the positive-frequency function obtained can be written in a more compact 



form with the ze-prescription The unfolded form presented in ( 71 ) allows us to see better 



the structure of its singularities. It is clear that all the other Green functions can be obtained from 
the positive-frequency function. 

Such an expression for the main contribution to the positive-frequency function has several 
physical implications. First, we see that the wave packet of scalar particles of the mass m in a 



slowly varying gravitational field behaves like a massive particle with the mass rh given in (67). In 
particular, if we neglect the term in rh? proportional to a 2 and consider the vacuum solution to 
the Einstein equations then the mass squared acquires the shift which is equal to — /i 2 /4 (for an 
analogous but not the same effect in quantum electrodynamics see |14|). This is a tiny negative 
quantity and is of the order of the Unruh temperature squared. One could expect the appearance 
of such a correction to the mass already from the expression ^ of the Fourier transformed Klein- 
Gordon operator. Notice that such a correction to the mass squared is also necessary to reproduce 



the standard asymptotic expansion of the heat kernel (18) and the propagator (see Eq. (96) below). 



This correction makes no trouble for the massive particles, but for the massless scalar particles it 
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seems result in the tachyonic dispersion law in the adapted system of coordinates and to the 
instability of a vacuum. If m 2 < and small, the additional correction appears in the positive- 
frequency function 



SDM(x,y)\ 



(c) 



-AV2( X) y) 



m 



(72) 



which comes from the accurate evaluation of the integral ( 70 ) and should be added to the expansion 



in the third line of (71). Certainly, we should not regard the expression (71) too seriously for the 



modes with the energies squared of the order of h 2 where the tachyonic effects may appear. These 



modes are out of the range of the applicability of the approximation we made in deriving (71). 



However, for the massive particles and for the high energy modes of the massless ones satisfying 



(50), (68) the positive-frequency function in the representations (70) or (71) with the corrections 



( 62 ) and ( 65 ) describes a propagation quite well. 

Second, the effective interval X 2 does not coincide with the interval (20). Indeed, it is easy 



to see for the worldline of a particle staying at rest, a = 0, that the equation X 2 {t) = has the 
nontrivial complex solutions (see Eqs. (20), (|24j)), while the equation X 2 {t) = has not. If we 
assume 



|r 2 | > l^(2/ 2 + £ (0) )<t/24|, 



then 



X 2 « ft 2 + ^HV 2 + E^)a + a 2 . 



(73) 



(74) 



The latter expression differs from (20) only by the terms in the square brackets which contain T. 
It readily follows from this observation the expansion of X 2 on an arbitrary worldline: 



X 2 
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4r 4 + ..., 



(y u u) 2 -[{al m + a^) 2 -a 2 l 



(75) 



where a^ m := (fcu) f uv'uF' and aff := (£u)(gu)hu are the gravimagnetic and gravielectric forces 
divided by the mass of a particle, respectively (see Appendix |A| for details). In particular, the ac- 
celeration squared afj is zero for a rest particle when u^ 1 = £^|£ 2 | -1 / 2 . We shall call the acceleration 
ajj as the Unruh acceleration since it determines the Unruh effect for a detector moving in a curved 
spacetime. 

As well known (see, e.g., |H [23]), the transition rate of the detector to the excited state is 
determined by the Fourier transform of the positive-frequency function taken on the worldline. 
More precisely, for the two points x(t\) and xfo) we introduce f := [t\ + T%)j2 and r := t\ — Tz, 
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and write £)( + )(x(ti), x(t 2 )) = D^ + \t,t). Then the rate of transition of the detector to the excited 
state d > is proportional to 



/ 



dTe- tu > T D (+ \f,T), (76) 



where the integration contour in the r-plane passes slightly below the real axis. The positive- 
frequency function D^ + \t) is singular at r = 0, where it possesses the pole of second order (for the 
non-isotropic worldline). All its other singularities are arranged symmetrically with respect to the 
imaginary axis. Therefore, up to the fourth order in r, they lie only on the real or imaginary axes 
when afj < and afj > 0, respectively. In the former case, the Unruh detector feels nothing since 
we can close the integration contour in the lower half-plane. Of course, there may be singularities 



that are not catched by the expansion (|75|), but they lie far from the real axis and so they are 



suppressed. In the case a? 7 > 0, the Unruh detector will detect the excitations at the reciprocal 



temperature 

(to = — , (77) 
au 

where we have replaced the factor \/l2 by 2ir matching the formula for the Unruh temperature 
with its flat spacetime analog for the hyperbolic motion. A concrete value of this factor depends 
on a character of the particle motion and requires a more detailed information about its worldline. 



Nevertheless, formula (77) provides a good approximation for the temperature by the order of 



magnitude. Thus we see that the Unruh detector can be employed to test the variations of the 



field Notice that if we used the interval (20) in order to describe the Unruh effect, we would 
obtain (V u ii) 2 for the Unruh acceleration and, consequently, arrive at the unphysical result that 
the resting Unruh detector gets excited in a stationary gravitational field. We should stress the 
important difference between the Hawking [3j and the Fulling-Unruh [H E] effects. The Hawking 
particle production is caused by the action of the gravitational forces, and the detector at rest may 
records this process. Contrarily, the Unruh detector responds to the action of the inertial forces. 
The latter can be defined as a 4-vector (see Appendix [Aj as long as the vector field exists on 
the manifold. 



The fact that the effective interval ( 74 ) differs from the geodetic interval results also in a nonva- 
nishing local expression for the radiation reaction force acting on the charged particle in a free fall. 
It is not difficult to show (see, e.g., |5-8j) that the retarded Green function constructed from the 



propagator with the interval ( 20 ) (see Eq. ( 96 ) below) yields the local contribution to the radiation 
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reaction force proportional to 

V u V u n M - (V u m) 2 m^, (78) 

on the vacuum solutions to the Einstein equations. This contribution vanishes for the geodesic 
motion and is not zero for a particle at rest. As for the Green function with the effective interval 



(75), the situation will be reversed: the radiation reaction force acting on the particle at rest will 
be zero, and it will be nonzero for the particle in a free fall. In a certain sense, the additional terms 
appearing in the effective equations of motion of a particle and stemming from the difference of 



( 71 ) from the DeWitt anzats take effectively into account the so-called tail term of the radiation 
reaction force. We postpone a thorough investigation of this problem to a future research, but 
should note that the charged particles can be also used to detect variations of the vector field 



One may wonder why the retarded (advanced) Green function following from ( 71 ) does not 
satisfy the general theorems concerning the structure of singularities of the fundamental solution 
to the hyperbolic partial differential equation |35[ I36j . The singularities of such a solution must 
lie on the characteristic cone determined by the geodetic interval, but we saw that the effective 



interval (74) does not coincide with this interval. The answer is that the expression for the main 



contribution to the positive-frequency function (71) is valid only for the point separation much 



larger than the wavelengths of the modes in the wave packet (see Eq. (68)). And so, the general 



theorem does not apply to the expression ( 71 ) . In order that the wave packet propagates along 
the characteristic surface, it should be delta-shaped in the space at the initial moment and consists 
of the ultrarelativistic modes. However, the delta-shaped wave packet is infinitely broad in the 



frequency space and does not obey the restriction (68). From the uncertainty relation for massless 
particles 

AujAx > 2vr, (79) 
where Aw is a width of the wave packet in the frequency space and Ax is its width in the space, 



we find that the above expression ( 71 ) for the Green function holds only for the wave packets with 



27T 

Ax»— , (80) 
where uq is the central frequency of the wave packet of the ultrarelativistic (massless) particles. 



Roughly speaking, the deviance of the singular surface of the contribution ( 71 ) to the positive- 
frequency function from the geodetic light cone just indicates that the wave packet with the ex- 
tension much larger than the wavelengths of its modes moves slightly different than an ideal point 
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massless particle. In the next subsection we shall see that in the limit (t, a) — > (violating the con- 



dition (68)) the structure of singularities of the positive-frequency function (52) is such as dictated 
by the theorems. 



B. The s-representation 

Hitherto we have analyzed the cj-representation of the positive-frequency function and derived, 
in particular, the expression for the major contribution to it using this representation. It turns 



out that under the assumption (50) we can take the integral over oj in (52) making the additional 



reasonable approximations. Thereby we shall derive the s-representation of the positive-frequency 
function. This representation will be employed to obtain the standard asymptotic expansion of the 
propagator in terms of the geodetic interval squared 2a at (t, <r) — > 0. 



If the condition (50 ) is satisfied, we can expand Sq staying in the exponent in formula (52 ) in the 
series in ui and retain only the terms at ui and oj . This can be easily done under the assumption 
that 



(81) 



The latter relations are valid in the weak field limit as well as for the spherically symmetric metrics. 
Then Sq in (52) keeps its form with the replacements E -4 E^ and b — > b^°\ while M 2 acquires 



the correction. The expression for this correction is rather huge. At the coincidence limit it reads 
as 



SM 



Vo = b<»E$b<» - lbME-]E {2) E-]bW + J|-(cos S xch s/ - 1)£-}&C°)- 



1 



b^E7^ + b^E^E m (E^f + i)^ 



Ctg SK 



ch sf 



sin sk 



)E^b(°\ (82) 



where x is defined as before, but with E — > E^°\ At the small s, this correction behaves as 

aE^a 



5M Z 
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+ 0{sa). 



(83) 



This expression is valid even if Eqs. (81) do not hold. The preexponential factor in (52) does not 



change apart from the redefinition of x. Then the integral over uj is reduced to the Hankel function 

,2, 



o \ a J i 



(84) 



where m 2 := M 2 + 5M 2 , the square root in the argument of the Hankel function has the cut along 
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Figure 1. The typical contour plot of the imaginary part of the function staying in the argument of the Hankel 



function in Eq. (87 1. The imaginary part changes its sign when reflected in the real axis. The insets depict the 
structure of singularities of the types (a) and (b) and the cuts near them. The thin lines on all these plots are the 
lines of the steepest descent. On the left panel: The thick line going along the real axis is the initial integration 
contour. On the right panel: The structure of singularities near the origin is depicted. The dashed line shows the 
steepest descent contour. The ends of the cuts lying on the real axis are the zeroes of a(s) nearest to the origin. The 
small line connecting two cuts intersects the integration contour at the saddle point. 



the positive real semi-axis, and a := s(t- — So), 
a = s 2 (C 2 + \bE~ l b) - s{t_ + l -bE- l fa)- 



- a— ctg sxo + {-- bE- 1 ) — (ctg sx - -. ) (- + E b). (85) 

4 A L Sill S% A 

In our case the Hankel function is expressed in terms of the elementary functions 



+ 2fha 1/2 



2ima 1 / 2 



a J d/2\— ) 2vr 1 /2 a 3/2 
Thus the positive-frequency function can be written as 



(86) 



(87) 



2|e 2 | 1 /4(x)(4 7 r) d / 2 |e 2 | 1 / 4 (y) 

The singularities of the integrand are located at the points (a) and (b) discussed above and also 
at the points where a(s) = 0. The cuts stemming from the square roots in the Hankel function 
and the determinant stretch between these singular points as depicted on Fig. [T] Recall that the 
integration contour in the s-plane lies a little bit lower than the real axis. 
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Unfortunately, the integral (87) cannot be evaluated exactly. Therefore we shall calculate it 
by the steepest descent method. To this end, we have to deform the integration contour and 
rise it from the forth quadrant to the first one (see Fig. [TJ the right panel). This gives the 
additional contributions to the integral from the cuts located on the positive real semi- axis, which 
are responsible for the vacuum polarization effects. We have already evaluated such contributions 
in the previous subsection and so here we concentrate on the major contribution to the positive- 
frequency function and its singularities. 



The main contribution to the integral ( 87 ) comes from the extremum positioned near the point 
s = 0. The singularities of the positive-frequency function appear when the two branching points 
nearest to the origin (zeroes of a(s)) approach the saddle point and pinch the integration contour. 
At (i, a) —> 0, this saddle point and the two branching points tend to the origin. In this limit the 
correction to M 2 is small and fh can be regarded as a constant. Hence, these saddle and branching 
points are determined solely by the function a(s). In the vicinity of the extremum we introduce a 
new variable z such that 

z 2 + z 2 = -a(s), (88) 
where zq is a constant to be determined and z = corresponds to the extremum point. Making 



use of the expansion (66), we have 



(J 9 ^ 1 „ _ S r, 9 1 , _o _ t-^{C\\\'? _ 1 



a = "T 



+ st- s 2 E - *-Wfa + b - [b 2 Q) + ^a(f + lE^fo - ^E 2 Q) a} + 0(s% (89) 



where we imply that t has a small negative imaginary part. Now the extremum so of the function 
a(s) can be found perturbatively assuming that t ~ a ~ / — > 0, i.e., we suppose that / is the smallest 
scale in the problem. Simple calculations give 

t 3 est 2 / 9 Co C4 s o 

2E 8 £3 V 16 £ 5 4E4 J > V ) 

where C3 and C4 are the coefficient of the expansion of the function a(s) at s 3 and s 4 , respectively. 
At the saddle point the function a(s) is proportional to the geodetic interval squared 

«(so) = -\{a 2 +ft 2 + ^H 2 [a{f 2 + C\V i h J -h i h j ))a-2hfat-h 2 t 2 }} (91) 

This determines the constant zq. 



Carrying out the calculations, we shall keep the terms in (87) up to the second order in I only. 



This will allow us to find the singular and finite parts of the positive-frequency function at I — ¥ 0. 
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In particular, according to this approximation m should be replaced by M as seen from (83) and 



(|87|. Thus we have 

B o ~a 2 ■■ c 3 +4c 4 s ... 3c 4 3 2 
E - 3c 3 s - 6c 4 s5, 2f» ; , z« ; —73(03 + 4c 4 s ) , 

z _ z (92) 

z c 3 + 4c 4 s 2 1 r 1 ^ / 1 /) \2i ^ 
s « so + t + ^ z + [c 4 + — (c 3 + 4c 4 s ) J 7775, 

where the overdots denote the derivatives with respect to s taken at s = sq. Introducing the 
notation Aq, A±, and A2 for the coefficients of the expansion in s of the preexponential factor in 



(87), we can write 



ds 2 A + Ais + A 2 s 2 Q 

— (^0 + MS + A 2 S ) « ; ~ + 

az z 

r . . , . , , C3 + 4caSq Z , . 3ca . . Z 2 . , 

+ [A x + 2A 2 s + (A + Ais ) - 2 — ] T2 + (A 2 + ^Aq)^. (93) 
The terms at the odd powers of z can be omitted as long as the integration over z is carried out on 



the symmetric interval. As a result, the integral in (87) can be cast into the form 



d [°° l A + A 1 s + A 2 s 2 3c 4 ,z 2 

ciz : + (A 2 + — Job 



-2Myfx 



2 , y 2 



2 _L ^2 



Z^ + Z, 







— (A + A lSo + A 2 s 2 )Ki(2Mz ) + ^(A 2 + ^A )Ko(2Mz ). (94) 
zzq z j 2z z 



Expanding A 1 / 2 and £ 2 in (87) near the midpoint up to the second order in I, we arrive at 

D {+) {x,y) « _^{[l+l J R^(^+^)(tr+^)]— ^i(2Mz )-[f / 2 +/ l 2 -?V A / l A ]ET (2Mz )}. 

327r z I 12 zo 6 3 J 

(95) 

Finally, we expand the Macdonald functions at the small argument |37j 

where 7 is the Euler constant and, as before, a is the world function, is its derivative. 

We see that the singularities of the retarded and advance Green functions lie on the geodetic 
light cone as it should be according to the general theorems. The expansion of these functions in 
terms of a coincides with the standard one (see [21, 23j) and is independent of Whereas the 
imaginary part of the propagator depends on the Killing vector nontrivially and "remembers" the 
vacuum state with respect to which the propagator is defined. We also see that the correction to 



the mass squared is reproduced again. Without it, in particular, the factor at the logarithm in (96) 



would not have the the covariant form independent of The mention should be made that the 



positive-frequency function (96) cannot be employed to analyze the Unruh effect or the radiation 
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reaction problem. For these problems the main contribution to the corresponding integrals comes 
from the sufficiently large point separation (smaller than L, but larger than the typical wavelength) 



where the expansions bringing us to (96) are not valid. The latter expansions may be justified in 



analyzing these problems only in the case of the extreme accelerations a caused by the external 
(nongravitational) force. In this case, the relevant contributions come from the point separations 
of the order of I ~ a -1 and the other terms entering the Unruh acceleration (75) can be neglected. 



Concluding this section, we write out the s-representation for the one-loop correction to the 



effective action induced by one bosonic mode. Comparing the formulas (52) and (60), one readily 
deduces 



^(1) 
lb 



f dx^3 t 2 Z)(+)(i,x,x) 



-iA-i) 



(97) 



where A is the ultraviolet energy cutoff. From the integral representation of the Hankel function 



(84) we observe that the differentiation with respect to t just results in the change of the overall 



sign and the replacement d — > d + 4 in (87) apart from the factor Atris. Consequently, we have 

i d+l dss 2 



r 
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dxy \g\ 



2 |£2|3/2( 47r )d/2 
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sin sk 
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1/2N 



(98) 



where mo equals m at a = and a(s) is given by the expression (85) with a = and t 



-zA" 



V. DISCUSSION 



In the previous sections we obtained an approximate solution to a purely mathematical problem 
of finding the Green function for a scalar quantum field on a stationary background under certain 
assumptions on the external fields (the Boulware vacuum is implied). We have shown that the 
two-point Green function of a scalar field on a curved background and the induced effective action 
(in particular, its imaginary part) depend nontrivially on the Killing vector i.e., this Killing 
vector field enters explicitly into the expressions and these expressions cannot be rewritten in a 
local form in terms of the 4-metric and its curvature alone. This dependence stems from the form 
of the Schrodinger equation for the quantum fields on a curved background 



d^\/\g\ 



1 



(99) 



where ^ is the state of quantum scalar fields, t is the time coordinate, and is the (regularized) 
operator of the energy-momentum tensor in the Schrodinger representation. The second term in 



the right-hand side of Eq. ( 99 ) is needed in order to throw the time derivatives to one of the fields 
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(j) entering the energy-momentum tensor. Then the Hamiltonian is diagonalized in the basis of 
the mode functions associated with The Schrodinger equation explicitly depends on the field 
g^ and the stationary states of this equation, in particular the vacuum state, also depend on it. 
For a flat spacetime with the vacuum state corresponding to the zero eigenvalue this dependence 
disappears, at least formally, because of the global Poincare symmetry of the Schrodinger equation. 

Assuming that the energy-momentum tensor of the matter fields is covariantly divergenceless in 
the case of a nonstationary metric |38| I39j. we arrive at the dynamical equations of a hydrodynamic 
form [9j for the vector field £ M with appropriate boundary and initial conditions. The form of these 



equations strongly suggests that the field g„ should be quantized. In this case, Eq. (99) is explicitly 
generally covariant as easy to see integrating it over t. The quantization of the hydrodynamic 
equations is more or less known [40j and is used in describing the cosmological fluctuations (see, 
e.g., [41, 42J). The quantum field g^ is decomposed into the condensate part and the fluctuations 
Sg^. Then the latter is represented in terms of the creation and annihilation operators, the particles 
being identified with the phonons. The properties of these phonons follow from the fluid energy- 
momentum tensor. In particular, as found in [S], the sound speed squared of such phonons equals 
—ifN/ot in the weak field limit, where <px is the Newtonian potential and a is the power of decay of 
the induced energy-momentum tensor at large distances from a gravitating body, i.e., T^ v ~ r~ a . 
In terms of the Feynman diagrams, the phonons 5g^ appear in the graphs. To say figuratively, it 
looks like the particles ring in moving in the condensate of the field g^. 

Although we do not touch in this paper the problem of quantization of the gravitational field, 
it should be noted that one or another solution to the time problem is necessary for the gravity 
quantization programme (for a review, see |43|). The vector field £^ and its dual, when integrable, 
can be employed to define the so-called world time [23] and provide the possible preferred definition 
of time. Such a method to solve this problem is not new and is known as the reference fluid 
approach (see, e.g., |43H46| ) . A distinctive feature of the approach proposed in [9] is that we do 
not introduce such a fluid by hand, but deduce its equations of motion from the requirement of 
the covariant divergenceless of the energy-momentum tensor (for other similar approaches see, for 
example, j37~50j). One of the purposes of the present paper was to show that the vector field 
defining such a fluid is already contained in the formalism of quantum field theory on a curved 
background. The existence of this vector (dynamical or not) allows one to define formally the 



gravitational and inertial forces as the 4-vectors (see |24[ 130] and Eqs. (A5), (A6) below). Notice 
that for certain Euclidian approaches (see, e.g., |51H53j ). where the Killing vector does not appear, 
the induced effective action contains the terms of the form RV~ 2 R which are essentially nonlocal. 
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These terms involve the inverse powers of the covariant d'Alambertian and the correct definition of 
such constructions requires the assignment of the initial and boundary conditions. In fact, this is 
equivalent to the introduction of the new fields into the model which also contribute to the covariant 
divergence of the energy-momentum tensor. 



Observe also a remarkable feature of Eq. (99) - it is invariant with respect to dilatations of the 
Killing vector. This demonstrates the invariance of a theory with respect to a possible choice of the 
energy unit. If we introduce the world time associated with the Killing vector (when it is possible) 
then this variable should be also stretched in such a way that the product "energy" x "time" will be 
invariant. Therefore we shall call this symmetry as the energy-time symmetry. One may check that 
all the observables appearing in this paper are invariant under such a transform. However, if the 
theory possesses an intrinsic scale, say, the cutoff scale, this invariance is violated. Formally, the 



violation follows from the fact that the right-hand side of (99) must be regularized. As a result, 
the additional terms involving appear and break the symmetry. The typical terms violating 
the energy-time symmetry are presented in Eq. (21) of j9] (see also |54|). One can distinguish 
two types of such violating terms: hard violating and soft violating. The former scale as certain 
powers under dilations of the Killing vector, while the latter do it logarithmically. The coefficients 
at the soft violating terms are well known and nothing but the scaling functions (/3 or 7) of the 
renormalization group taken at that scale where the model is in a perturbative regime (if it exists). 
Note that all the particles of the model, even with large masses, contribute to this coefficients on 
equal footing (for a concrete calculations in the standard model see |55]). The relation of these 
coefficients with the scaling functions is not surprising (see, e.g., |10[ lll| ) as long as the energy- 
time symmetry involves a dilatation. Hence, if the theory were conformal on a quantum level, 
such anomalous terms would not arise. The restriction on the parameters of the models implying 
/3 = can be found, for example, in |56j . According to [D], the anomaly of the effective action 
under the energy-time dilatations defines the enthalpy density of the condensate g^, although the 
hydrodynamic equations for g^ are not empty even in the case of a vanishing enthalpy. 



Appendix A: (3 + 1) with the Killing vector 

In this appendix, we collect some formulas regarding the differential calculus on the Riemannian 
manifold with the Killing vector. If the metric possesses the Killing vector £^ then the following 
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useful relations hold 



fag X = 0, gxh x = 0, g x Vxfa = 9{f.f*]\h X , 
= g 2 h 2 , V x fag<+ = -\f\ V u g u = \fa- h [p g u) , 



(Al) 



and for the curvature 

1_ . 1 



R^g v = \f - g 2 V x h x , g x R Xu = f u \h x - ^V x fa - g„V x h x , 



where fa = d^g u ] and = d n In . 

Using the Killing vector we can construct the projected connection, 



(A2) 



O := r«j + f {f a {il m + 2h a gii g p - 2g a h {m ), (A3) 



such that 



(A4) 



V^p ff = 0, = 0, V n g v = - fa, 

V»9 P a = y/ M p5<x) + 2£ 2 / W ., V ^ = -^-f^g a \ V p g^ = 0. 

Then the equation for geodesies of the metric g uv takes the form (cf. [24J) 



V u u» = V u u» - (tu)(r u u» + (guW - 2{hu)g») = => 

g^ u u u = V u (g nu u u ) - {iu)(f pv u u + {gu)h^) = 0, (A5) 

The second term on the left-hand side of the first equation can be regarded as the acceleration 
caused by a stationary gravitational field [30J. The first term in parentheses is associated with the 
so-called gravimagnetic force, whereas the second term is responsible for the gravielectric component 
of the force. The factor at these parentheses is the energy of the particle divided by its mass. In 
particular, the vector field £^ allows one to define formally the inertial forces as the 4- vector 

fLr ■= + (WX + M^) " 2{hu)g»] & tf ner + f£ + f» = 0, (A6) 

where f 1 * := m\7 u u^ and fg := m(^u){f^ u u u + {gu)h^ — 2{hu)g^) are the external (nongravitational) 
and the gravitational forces, respectively. 
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The Riemann and Ricci curvatures of the projected connection read as (see also 

+ i 2 {g[ah^h v g p - g[ a h^h p g v + ^h [a f p]u g p - ^g [a f p]u h p + g [a V ^h v g p - g[ a V p] h p g u - 

- U(Kg p - h p g v ) - f u p h w9fl] - \h [(T f p g v + \g { J^h v ) + ^<&,<7 [ff /„f (// + 2h x g p ), 

R„u = R„u + KK + V„hu - %2~\U P fv + 9{»S P f v P + 2 9ll g u (h 2 - V p h p ) + g {p f v p h p ] + jg^f, 

(A7) 

where f 2 := f^f^. Such a notation is used for scalars. As for matrices, f 2 := fu,xf\, and then 
Sp/ 2 = — f 2 . The scalar curvature becomes 

R = r v R„u = R + 2V p h p + ^f 2 . (A8) 

Besides, we have 

V% = r u V tl h u = V x h x -h 2 , h { h l = h 2 , fijf J = f 2 . (A9) 
Also we need the quantities in the adapted coordinates in the weak field limit |24] 

Rikim ~ (2<%£ m ]fc — 3njn[;(5 m ] fc + 3n fc n^ m ]j), R km m — ^ (S km — 3n k n m ), 

^2 ^.2 

RikRj ~ + SniUj), R ik i m Rj m w (55y — 3^^), 

3r 2 r 3r 2 3r 2 

V 2 R lj dx l dx j = (1 - -fY 1 ^ 2 - ^(r 2 d9 2 + r 2 sin 2 9dip 2 ) « -^-§(<% - 3mn j )da^dx i , 

-'' ' '" ''" (A10) 

where nj = x«/r and V-Rjj is calculated for the Schwarzschild metric. Recall that = and 
V x h\ = for this metric. 



Appendix B: Some expansions 

In this appendix we give the expansions of some quantities appearing in the main text. For the 
metric in the Riemann normal coordinates we have |27| [29] 

9ij = Sij - -R ikj iy k y l - ^ s RikjW k y l y s + ^{—R^Rukm - 6V ' S nRikji)y k y l y s y n + 0{y 5 ), 
-g = i- \Ri j y i y j - ^ l R ]k y i y j y k - (^i^ + ^R m%m R m k i n - \Ri 3 Rki)y i yh k y l + 0( y 5 ). 

(Bl) 



Also we deduce 

ie 2 r 1/4 (x)ie 2 r 1/4 (y) = ie 2 r 1/2 (p)(i - + o^), 

i 

6 J 



A-^x.y) = (1 - -RijG^ + O^)). 



(B2) 
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In order to prove the relation ( 24 ) , the following expansions may be useful 



-1 = eV) 2 + UiV? + 2[£ 2 ((u ) 2 h i u i + u°u°) + u i « i ]r+ 

+ [t 2 [{u®) 2 {hiij + 2(/j^) 2 + Vihjtful) + 4u°u°/i i u i + (u ) 2 

+ y frjuW] + (u') 2 + u^Jr 2 + 0(r 3 ), 
(V u n) 2 = e 2 (^°) 2 + + £ 2 [eV) 2 (/> V) 2 + 2u°h%v? - u%v?) 
- 2u°(u l f ijU j + vPhitf) + AhitfvPu + 4(/i^) 2 (u ) 2 ] + O(r), 
C 2 (P) = e 2 [l + MV + [hiij + (/i^ 4 ) 2 + Ivi^-«V)y] + 0(r 3 ). 



(B3) 



Here all the fields are taken at the origin of the system of coordinates described after Eq. (21). 
Notice that u° = g^u^ in this frame. 
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